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Talk Outline 


1. Ultra-spectral IR Instruments 

2. IR-only Retrieval Algorithm Introduction 

3. Retrieval Simulation Analysis 

4. Retrieval Demonstration from Real Satellite Measurements 

5. Validation and Inter-comparison (JAIVEx) 

6. Future Ultra-spectral IR Instruments on GEO Satellites 

7. Summary 
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NAST-I and S-HIS: Since July 1998 



NPOESS Airborne 
Sounder Testbed 
(NAST) 

on ER-2, WB-57, 
Proteus aircraft since 
July 1998 



iR-2 Aircraft 


NPOESS Airborne 
Sounder Testbed 
Interferometer (NAST-i) 


Electronics, Control, 
faavigation, Processing E 
Recording 


NPOESS Airborne 
Sounder Testbed 
Mkrow$v% Red&omotor 
Sounder (MAST-M) 
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AIRS: Since May 2002 


Atmospheric 
InfraRed Sounder 
(AIRS) instrument 
(by NASA) on Aqua 
Satellite launched on 
4 May 2002 
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IASI: Since Oct 2006 


Infrared Atmospheric 
Sounding 

Interferometer (IASI) 
instrument (by CNES 
/ EUMETSAT) on 
MetOp-A Satellite 
launched on 19 
October 2006 
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Retrieval Parameters from These Systems 


Brightness Temperature or Radiance Spectrum 
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Retrievals under clear conditions: 

• Surface properties (skin temp and 
emissivity). 

• Atmospheric temperature and moisture 
profiles. 

• Atmospheric CO and 0 3 abundances. 


Retrievals under cloudy conditions: 

• Atmospheric profile through optically thin 
cirrus clouds and above optically thick 
clouds. 

• Effective cloud parameters (i.e., cloud top 
pressure, particle size, and optical depth). 
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LaRC IR Retrieval Algorithm 


PART A: REGRESSION RETRIEVAL (Zhou et al., GRL 2005) 

Using an all-seasonal-globally representative training database to diagnose 0-2 cloud layers from 
training relative humidity profile: 

A single cloud layer is inserted into the input training profile. Approximate lower level cloud using 
opaque cloud representation. 

Use parameterization of balloon and aircraft cloud microphysical data base to specify cloud effective 
particle diameter and cloud optical depth: 

Different cloud microphysical properties are simulated for same training profile using random 
number generator to specify visible cloud optical depth within a reasonable range. Different habitats 
can be specified (Hexagonal columns assumed here). 

Use LBLRTM/DISORT “lookup table” to specify cloud radiative properties: 

Spectral transmittance and reflectance for ice and liquid clouds interpolated from multi-dimensional 
look-up table based on DISORT multiple scattering calculations. 

Compute EOFs and Regressions from clear, cloudy, and mixed radiance data base: 

Regress cloud, surface properties & atmospheric profile parameters against radiance EOFs. 

PART B: 1-D VAR. PHYSICAL RETRIEVAL (Zhou et al., JAS 2007) 

A one-dimensional (1-d) variational solution with the regularization algorithm (i.e., the minimum 
information method) is chosen for physical retrieval methodology which uses the regression 
solution as the initial guess. 

Cloud optical/microphysical parameters, namely effective particle diameter and visible optical 
thickness, are further refined with the radiances observed within the 10.4 pm to 12.5 pm window 
region. 
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LaRC Algorithm Flowchart 



Hybrid Retrieval Algorithm Flowchart 



Calibrated 

spectral 

radiances 


EOFs & regression coefficients 
with categories: (1) clear, (2) all 
clouds, (3) ice cloud grouped with 
cloud heights, and (4) water cloud 


Regression retrieval using “all 
cloud” coefficients category to 
predetermine cloud height 



Historical training data 
with cloud parameters for 
radiance simulation including 
a realistic cloud radiative 
transfer calculation 


Regression retrieval using “cloud 
height grouped” coefficients for ice 
to iterate and finalize the retrieval 


Regression retrieval using “cloud 
height grouped” coefficients for 
water to iterate and finalize the 
retrieval 



Regression retrieval using “clear” 
coefficients 


Cloud parameter fixed for 
refining thermal properties. 


Cloud optical thic kn ess and particle size 
are refined with cloud effective micro- 
window radiances: 11.1 pm to retrieve 
cloud optical depth, and 11.1 pm & 12.5 
pm windows to retrieve cloud particle 



Compute Jacobian Matrix: analytical 
Jacobian for thermodynamic parameters; 
and numerical Jacobian for cloud 
parameters (i.e., x and P cld)- 

Matrix Inversion Solution: the 
regularization algorithm (or minimum 
information method) in a nonlinear least- 
squares spectral radiance fitting mode. 


5R = A5Q, where A- = oR/oQj 

SQ n = [(A^EA+yiy^E-\5R+A8Q)] nA 

5R n = R 1T1 -R c (Q n ) iterate to minimum. 

5Qn = Q n -Qo (but 5 t = O.lx, 5P cld = P cld (k)-P cld (k-1)] 


n = iteration number; k = pressure (P) grid number. 

Q = variables [i.e., T(P), q(P), Ts, Ps, Pcld, x, De(x) ...] 
Q 0 = regression retrieval. 

Rm = measured radiance 

R C (Q) = calculated radiance from Q 

A = Jacobian Matrix 

Ay = [R i (Q J )-R i (Q J -5Q J )]/ 5Qj (for cloud parameters) 
E = error covariance matrix 

y = Lagrangian multiplier (optimized in the iteration) 






Statistical EOF Physical Regressions 


Simultaneous 1-D Var. Iterative Matrix Inversion 
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IASI Clear Retrieval Analysis 
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Synthetic analysis: the truth profile (i.e., the radiosonde observation) is known and the 
retrieval can be directly compared with the truth to define retrieval accuracy due to (1) 
instrumental noise and (2) retrieval error introduced mainly by so-called “ill-posed” 
retrieval model. The disadvantage of this approach is that forward radiative transfer 
model error is not included. 
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IASI Cloudy Retrieval Analysis 
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T Variance of Test Dataset and IASI Retrievals 


Under Clear 
Conditions 


Under Cloudy 
Conditions 
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IASI Retrieval Demo: Cloud Parameters & Ts 
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IASI vs. GOES-12: Cloud 
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RH Horizontal Distribution (%) 24.49 km 
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IASI Retrieval Demo: Moisture Distribution 
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Joint Airborne IASI Validation Exp. (JAIVEx) 



Location/dates: 

Ellington Field (EFD), Houston, TX, 14 Apr - 4 May, 2007. 
Aircraft: 

NASA WB-57 (NAST-I, NAST-M, S-HIS); 

UK FAAM BAel46-301 (ARIES, MARSS, SWS; dropsondes; in-situ 
cloud phys. & trace species; etc.). 

Satellites: 

Metop (IASI, AMSU, MHS, AVHRR, HIRS). 

A-train (Aqua AIRS, AMSU, HSB, MODIS; Aura TES; CloudSat; 
and Calipso). 

Ground-sites: 

DOE ARM CART ground site (radiosondes, lidar, etc.) 
Participants: 

include NASA, UW, MIT, IPO, NO A A, UKMO, EUMETSAT, 
ECMWF, ... 
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JAIVEx Case Study - Validation (2007.04.29) 




GOES-12 IR image 
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Retrieval Consistency Check: fitting residual 


Fitting Residual: STD of the difference between measured and retrieval 

simulated brightness temperature over physical retrieval 
channels. 


Clear Fitting Residual 




Cloudy Fitting Residual 
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i 4 Retrieval Consistency Check: Fitting Statistics 


Clear fitting sample (35.36N, 93.67W) 
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IASI Retrieval: Cloud Parameters 


GOES-12 IR image 
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IASI Retrieval: Surface Parameters 
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IASI Retrieval: ATemp and RH Fields 
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IASI Retrievals vs. Radiosondes 
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High-V ertically-Resolved Retrievals 
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IASI Regression vs. Physical Retrieval 



1. The retrieval improvement based on the EOF statistical regression through 
physical iterative retrieval is only contributed by IASI measurements as the 
minimum information methodology used. 


2. A high-vertically-resolved atmospheric structure is captured very well by IASI 
measurements and/or retrievals; not only in the troposphere, but also in the 
boundary layer. 
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IASI (15:48 UTC) vs. AIRS (19:30 UTC) 
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NAST-I: Connection between IASI and AIRS 
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Difference is mainly due to 

• instrument difference between IASI 
and AIRS, 

• spatial resolution difference between 
NAST-I and IASI (or AIRS), and 

• retrieval uncertainty including 
radiative transfer models difference. 
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Summary and Future Work 


1. A state-of-the-art IR-only retrieval algorithm has been developed with an all- 
seasonal globally representative EOF physical regression and followed by 1-D 
Var. physical iterative retrieval for IASI, AIRS, and NAST-I. 

2. The benefits of this retrieval are to produce atmospheric structure with a 
single FOV horizontal resolution (12 km for IASI and AIRS), accurate profiles 
above the cloud (at least) or down to the surface, surface parameters, and/or 
cloud microphysical parameters. 

3. Initial case validation indicates that surface, cloud, and atmospheric structure 
(include TBL) are well captured by IASI and AIRS measurements. Coincident 
dropsondes during the IASI and AIRS overpasses are used to validate 
atmospheric conditions, and accurate retrievals are obtained with an expected 
vertical resolution. 


4. JAIVEx has provided the data needed to validated retrieval algorithm and its 
products which allows us to assess the instrument ability and/or performance. 

5. Retrievals with global coverage are under investigation for detailed retrieval 
assessment. It is greatly desired that these products be used for testing the 
impact on Atmospheric Data Assimilation and/or Numerical Weather 
Prediction. 
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Fine-scale atmospheric horizontal features with high 
vertical resolution from satellite global observations 
with advanced ultra-spectral instruments have been 

realized for the first time. 
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k Geosynchronous Imaging FTS (GIFTS) 


Revolutionary Technology for Observing Atmospheric 
Temperature, Moisture, Winds, and the Transport of Pollutant Gases 


An Opportunity for Greatly Improved Environmental Forecasts 



Large area format focal 
plane detector arrays. 
Fourier transform 
spectrometer 
Geostationary satellite 
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GIFTS Science Objective Demo 


Wind Profile Objective: 

High resolution, vertical (1-2 km) and horizontal (50 km), cloud and water 

vapor wind profiles with a 4 m/s accuracy. 



An algorithm to derive 
clear-sky, altitude- 
resolved atmospheric 
motion vectors (AMV) 
is being developed and 
evaluated using 
simulated ultra-spectral 
data based on the fr 

L i 

GIFTS instrument 
[Velden et al., 2004; 
2005]. 
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Science Products: 

• Water vapor (soundings, fluxes, winds): 

• Temperature (sounding, stability): 

• Wind Velocity: 


e < 20% per 1 -2 km layers 
8 < 1°K per 1-2 km layers 
8 < 4 m/s per 2 km layers 


• Carbon monoxide concentration (2 Layers): 8 < 10% per 5 km layers 


• Clouds (altitude, optical / microphysical properties, “winds”) 

GIFTS Sampling Characteristics: 

• Two 128x 128 Infrared focal plane detector arrays with 4 km footprint size 

• A 512 x 512 Visible focal plane detector arrays with 1 km footprint size 

• Field of Regard 512 km x 512 km at satellite sub-point 

• 11 second full spectral resolution integration time per Field of Regard 

• ~ 80,000 Atmospheric Soundings every minute 


Ozone concentration (4 Layers): 
Surface Temperature: 

Aerosol Concentration and Depth: 


8 < 10% per 8 km layers 

8 < 0.3°K for sea; 8 < 1°K for land 

8<TBD 
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GIFTS Spectral Radiance Accuracy 


Spectral Radiance Accuracy Objective: 


• Spectral Coverage: 

• Spectral Resolution: 

• Spectral Stability: 

• Absolute Radiometric Accuracy: 

• Radiometric Noise: 
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* ; Components & Optical System Ray Trace 
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* ; GIFTS EDU Ground Based Experiment 
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Ground Based Experiment: OME 


To demonstrate GIFTS 4-D observation ability with the GIFTS EDU on the ground, 
an experiment called Outgas Monitoring Event (OME) is designed to monitor both 
temporal and spatial variations of the atmosphere with the GIFTS EDU when 
Sulfur Hexafluoride (SF 6 ) gas is released in the near field 
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Ground Based Experiment: OME 


To demonstrate GIFTS 4-D observation ability with the GIFTS EDU on the ground, 
an experiment called Outgas Monitoring Event (OME) is designed to monitor both 
temporal and spatial variations of the atmosphere with the GIFTS EDU when 
Sulfur Hexafluoride (SF 6 ) gas is released in the near field 
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Illustrations of SF6 outgas measurements provided by the GIFTS EDU on 18 September 
2006 - a demonstration of the GIFTS image tracking ability. 
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Ground Based Experiment: MTE 


One of the fundamental advantages of the GIFTS concept is illustrated by the 
spectral images shown in the video below, an example of the lunar and 
atmospheric images obtained with two GIFTS infrared detector arrays indicating 
that the lunar and atmosphere image is a function of wavenumber. 
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Illustrations of lunar and atmospheric measurements provided by the GIFTS EDU on 11 
September 2006 - a demonstration of the GIFTS imaging spectrometer measurement ability. 
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Fine-scale atmospheric horizontal features with high 
vertical resolution from satellite global observations 
with advanced ultra-spectral instruments have been 

realized for the first time. 

GIFTS EDU imaging and tracking capability 
demonstrates that quasi-continuous measurements of 
the moisture flux can be obtained for timely forecasts 

of storm intensity changes. 
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